Priedel phases and phases of transmission amplitudes in quantum scattering systems 
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We illustrate the relation between the scattering phase appearing in the Friedel sum rule and 
the phase of the transmission amplitude for quantum scatterers connected to two one-dimensional 
leads. Transmission zero points cause abrupt phase changes ±7r of the phase of the transmission 
amplitude. In contrast the Friedel phase is a continuous function of energy. We investigate these 
scattering phases for simple scattering problems and illustrate the behavior of these models by 
following the path of the transmission amplitude in the complex plane as a function of energy. We 
verify the Friedel sum rule for these models by direct calculation of the scattering phases and by 
direct calculation of the density of states. 



I. INTRODUCTION 
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The phase is an essential concept in quantum scat- 
tering theory. Some key results and techniques, suck 
as partial wave expansion,!!! and the Friedel sum rulea 
depend in an explicit manner on the scattering phase. 
The Friedel sum rule connects the density of states to 
the charge (or a charge difference) of the system ida the 
phase of the eigenvalues of the scattering matrix.Ba Since 
the Friedel phase is related to the density of states it 
is also connected to thermodynamic statistical mechan- 
ics quantitieaj'tj like the persistent current J3 In addition 
to the total density of states the scattering phases alsa 
play an important role in the partial density of statesQ 
(density of states with a preselection or postselection of 
the incident or exiting quantum .-channel) and in trans- 
port coefficients like capacitancesfl and charge relaxation 
resistances.!!!] 

The principal aim of this work is to investigate the 
behavior of the phase in simple scattering problems as 
they frequently occur in mesoscopic physics. In partic- 
ular we would like to understand which phases are con- 
tinuous functions of external parameters (Fermi energy, 
magnetic field or Aharonov-Bohm flux) and which phases 
arc permitted to exhibit jumps as a function of the exter- 
nal parameters. We consider coherent quantum scatter- 
ing systems connected to two semi-infinite, one-channel 
leads in the absence of a magnetic field and without spin- 
orbit scattering. Consider for instance the transmission 
amplitude t which determines the transmitted current 
amplitude if there is an incident current of unit ampli- 
tude. The transmission amplitude can be expressed in 
term of its modulus |i| and its phase 9 {t) , 
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If as a function of energy \t\ is always positive the path of 
the transmission amplitude will encircle the origin of the 
complex plane but always stay at a finite distance from 
it. An example of such a behavior is shown in Fig. nl As 
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FIG. 1. Transmission amplitude as a function of energy for 
a resonant double-barrier. 
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FIG. 2. Transmission amplitude as a function of energy for 
a wire with a side branch. 



a consequence the phase 9 {t) is a continuous function of 
the energy. If on the other hand |i| is zero for a certain 
value of the external parameter, the path of the trans- 
mission amplitude in the complex plain will necessarily 
pass through the origin. As a consequence the phase 9 {t) 
will jump at this energy by w. Such a behavior is shown 
in Fig. g. A specific model giving rise to the behavior 
shown in Fig. n is the transmission amplitude of a simple 
one-dimensional resonant tunnel barrier structure. The 
transmission probability never vanishes and the phase 
9 {t) is a continuous function of energy which increases 
in each complete revolution by 2ir. In contrast, Fig. 2 
shows the transmission amplitude for a simple model of 
a side branch of finite length attached to a perfect wire. 
The transmission amplitude as function of energy passes 
through zero and the phase which increases by 7r through 
each revolution also exhibits a jump of 7r (or — ir) bringing 
the phase back to its origin. 

The two models with the very distinctive behavior can 
be combined and the evolution of the transmission am- 
plitude of such a combined model is shown in Fig. 3. 
Now the graph shows very many revolutions through the 
origin and occasionally a revolution around the origin. 
We believe that this reflects the generic behavior of the 
transmission amplitude in the complex plane. 

Obviously the phase of the transmission amplitude in 
the second example, Fig. 2, since it is not a continuous 
function, can not play the role of the scattering phase 
which is used in the Friedel sum rule. For a system with 
a density of states p there must exist a phase which we 
denote by 9^> and which we call the Friedel phase, such 
that its energy derivative is directly related to the density 
of states 
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Since the density of states should be a continuous func- 
tion of the energy for the scattering problems we have in 
mind, the Friedel phase must also be a continuous func- 
tion of energy. One aim of our work is to investigate and 
illustrate the behavior of the different phases 9^ and 
0(f) and to investigate their connection to the scattering 
matrix of the problem. 

The problem investigated here is of interest in connec- 
tion, with experiments by Yacoby et al.lill and Schuster et 
al.E3 In the experiment of Yacoby et al. an Aharonov- 
Bohm ring with a quantum dot was investigated in a 
two terminal geometry. In a two terminal-geometry the 
Aharonov-Bohm effect exhibits a parityxQ As function 
of the Aharonov-Bohm flux the conductance is either a 
local minimum at zero flux (positive parity) or a local 
maximum (negative parity). In the experiment of Ya- 
coby et al. it was observed that over a sequence of more 
than a dozen Coulomb blockade peaks the parity changes 
at each peak in an identical manner. It is observed that 
the parity is positive to the left of the Coulomb blockade 
peak and negative to the right of the Coulomb blockade 



peak. Such a behavior is incompatible with a simple res- 
onant tunneling model Fig. pi but would be in accordance 
with the evolution of the transmission amplitude shown 
in Fig. 2. The second important experimental fact is that 
the phase drops by 7r between Coulomb blockade peaks 
as shown in the experiment of Schuster et al. Again this 
behavior is compatible with Fig. 2 but not with Fig. [y. A 
number of effortsll3c3 have been made to explain the fact 
that the behavior of the parity is the same at each peak. 
Ref. [l3] suggested a screening effect, Ref. [l5] alluded to 
degeneracies, Refs. [17lf2C| . |2l| and 24 proposed asymmetric 
deformation of the dot and repeated tunneling through 
exactly the same state. The same mechanism due to a 
dot whicb4s only semi-chaotic is supported by Silvestrov 
and Imrycia who point to the large variation in the lifetime 
of states in a semi-chaotic dot. A possible role of zero's 
of the transmission probability, Fig. 2, was suggested to 
us in informal communication by Levy YeyaAiE2Laad has 
since found interest in a number of works.li§E3S'E3~^Z 
We mention here especially the work by Ryu and ChoEi 
who investigate an AB ring with a dot which is also con- 
nected to a side branch. Apart from Ref. [13] the Friedel 
sum rule has-jthus far only found interest in the recent 
work of Lee.E3 The work presented here is in very close 
connection to the papers by Ryu and Cho and the paper 
by Lee. 
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FIG. 3. Transmission amplitude as a function of the par- 
ticle energy for a model which combines the resonant tunnel 
barrier and the side branch. 

It is possible that the experimental observation of iden- 
tical behavior over many subsequent conductance peaks 
is in fact generic (i.e. is observed for a fully chaotic quan- 
tum dot) and does require a system with special prop- 
erties which generates tunneling through the same state. 
The parity of the Ahronov-Bohm effect was observed in 
a multi-channel GaAs ring already almost a decade-ago 
in an experiment by Ford, Washburn and Fowlcr.E3 In 



this experiment the ring is coupled to a back gate and 
the experiments shows wide regions in the gate voltage - 
magnetic field plane in which the parity of the Aharonov- 
Bohm effect remains the same. 

In this work we investigate the various phases of the 
scattering matrix and their relation to the phases of 
the transmission amplitude and the Friedel phase 9 {t) 
comparing in detail two simple model systems, a one- 
dimensional resonant tunneling problem and a perfect 
wire with a side branch and briefly discuss also the com- 
bination of these two models. 



Apart form a constant tt/2 the Friedel phase is deter- 
mined only by the phase 9 of the scattering matrix. In 
particular, as also emphasized by Lee,E3 it would be in- 
correct to identify the Friedel phase with the argument 
of the transmission amplitude. 

The derivative of the Friedel phase 9 {n with respect to 
the particle energy E can be related to the energy deriva- 
tives of the scattering matrix. The density of states can 
also be expressed in terms of the scattering matrix.El This 
gives us a relation which connects the energy derivatives 
of the Friedel phase and the scattering matrix with the 
density of states, 



II. FRIEDEL PHASE AND TRANSMISSION 
AMPLITUDE PHASE 

We consider systems connected to two semi-infinite, 
one-dimensional leads. The system itself can be a multi- 
dimensional, i. e. a quantum dot. We consider phase- 
coherent scattering and neglect particle-particle interac- 
tions and inelastic scattering effects. In such systems 
the scattering matrix S in the particle energy E is rep- 
resented as a 2 x 2 matrix in which Sjj is the reflec- 
tion amplitude back into the j-th lead for carriers inci- 
dent in the j'-th lead (j = 1,2) and Sjk is the trans- 
mission amplitude from the k-th lead to the j'-th lead 
(j 7^ k, j = 1,2, k = 1,2). The scattering matrix is an 
unitary matrix, S^ 1 = S^ . This condition guarantees a 
conservation of the particle current. Furthermore it im- 
plies that the eigenvalues of the scattering matrix S are 
on the unit circle. Therefore the eigen value of the scat- 
tering matrix S can be represented as e 2 *^ with a real 
quantity £j (j — 1,2). Below we show that we obtain 
the Friedel sum rule if we define the Friedel phase 9 U) by 
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As for any definition of a phase, Eq. (ph, defines the 
Friedel phase only up to a multiple of n. A unitary 2x2 
matrix can be parameterized as follows, 



S = 
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with real phases 9, ipi, <p 2 and 0. The eigenvalues of this 
matrix Xj — exp(2z£j) are determined by 

sin(2£j — 9) = sin^i • cos0. (5) 

The eigenvalues are independent of ip 2 - Furthermore if 
£_ is a solution of Eq. (fsj) then a second solution is 
2£ + — 9 = it — (2£_ — 9). Consequently the sum of the 
two phases of these eigenvalues is 2£ + + 2£_ = 7r + 29. 
Thus according to Eq. (||) the Friedel phase is given by 
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Integration of this relation over the energy interval 
[Ex, E2] gives the Friedel sum rule, which thus states that 
the difference in phase 9 U) (E2) — 6 U) (E\) is equal to the 
number of particles N(E%, Ei) multiplied by 7r in the sys- 
tem in this energy interval, 



e^(E 2 )-6 if) (Ei)=7TN(E 2 ,Ei) 
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The Friedel phase is thus a a continuous function of en- 
ergy (either E 2 or E\) since the density of states is a 
continuous function of energy. 

Another phase which is frequently discussed is the ar- 
gument 9 (t) of the transmission amplitude defined by Eq. 
(|l|) . Below we will consider only systems with time rever- 
sal invariance. For such systems the scattering matrix is 
symmetric 6*12 = S21 and consequently ip 2 = or tp 2 = n. 
To be definite, we take, tp 2 — 0. In this case the phase of 
the transmission amplitude is given by 



7rO(cos< 



(9) 



where Q(x) is the step function of x. In contrast to the 
Friedel phase, the phase of the transmission amplitude 
is thus in general not a continuous function of energy, 
but exhibits a jump of it when cos <p changes sign. The 
Friedel phase and the transmission amplitude phase are 
not completely independent. From Eqs. (H) and (0) we 
find 
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Differentiating Eq. ( |l0|) with respect to energy (or any 
other variable of interest) we find, 



A9 (t) = A9 (f) ±7r£(cos0)AcoS( 
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where we have used the abbreviation A = d/dE. Eq. 
( Jll| ) shows that only if the condition cos (j) 7^ is satis- 
fied for all values of the parameters is the phase of the 
transmission amplitudes equal to the Friedel phase. The 



condition cos 0^0 means that the transmission prob- 
ability | t | 2 is nowhere equal to zero. Therefore it is 
the existence of zero points of the transmission probabil- 
ity which is at the origin of the difference of the Friedel 
phase and the phase of the transmission amplitude. 

To be more precise, in addition to a zero point in the 
transmission probability, it is also required that the en- 
ergy derivative of d cos <j)/dE is non-zero at the transmis- 
sion zero-points. To see this, we consider phase changes 
due to the energy E only, and assume that the number 
of zero points of the transmission probability is hnite or 
countable. The zero points of the transmission probabil- 
ity determines a sequence of energies which we denote by 
E {n \ n = 1, 2, ... Using Eqs. (|7|) and (III]) we obtain 
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where the function sgn(a;) of x is the sign function. It 
may be noted that the second term of the right-hand 
side of Eq. ( |l2| ) is zero even at a transmission zero point 
if d cos 4>/dE is zero at such a point. However, we expect 
that such cases are unlikely, and we proceed by assum- 
ing that d cos 4>/dE is not zero at any transmission zero 
point. 

To summarize: The Friedel phase is a continuous func- 
tion of energy. On the other hand, the argument of the 
transmission amplitude, .might exhibit jumps. These con- 
clusions agree with Lee.EZl 
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FIG. 4. Schematic representation of the path of the trans- 
mission amplitude as a function of the particle energy for the 
case that there is no transmission zero. The energy interval 
covered is that needed for the transition from one resonant 
peak (in the transmission probability) to the next resonant 
peak. The phase changes by about 7r. 



III. EVOLUTION OF THE PATH OF THE 

TRANSMISSION AMPLITUDE CONNECTING 

CONSECUTIVE RESONANCES 

We now discuss, from a general point of view, the con- 
ditions which lead to the behavior shown in Figs. EL g 
and a. We consider a set of parameters for which the 
transmission probability shows a series of peaks as func- 
tion of energy. We now ask: What portion of the path 
of the transmission amplitude shown in Fig. |l] and Fig. 
is traced out if we increase the energy from its value 
at a peak in the transmission probability to a value cor- 
responding to the subsequent peak in the transmission 
probability? Let us denote the energy of the n-th peak 
in the transmission probability by £ n . In addition, let us 
consider the condition 
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which according to the Friedel sum rule implies that one 
particle is added to the scattering region. Using this con- 
dition and integrating both sides of Eq. (O) with respect 
to the energy from £ n to £ n +\ we find 
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In case A there is no transmission zero point in the energy 
interval {£ n ,£ n +\), and in case B there is a transmission 
zero point in the energy interval (£„, £ n+ \). In case A the 
phase 9 {t) evolves by it through the consecutive resonant 
peaks, so that the resulting path of the transmission am- 
plitude is as shown in Fig. [|. On the other hand, in case 
B the phase 9 {t) of the transmission amplitude increases 
by 2ir (or 0) between the consecutive resonant peaks, so 
that the trace of the transmission amplitude is as shown 
in Fig. |. 

It is clear that Fig. 1 is composed of two paths of 
the type shown in Fig. 4, and Fig. g is the result of a 
path of the type shown in Fig. 0. Consequently, there 
is a profound difference in the behavior of these two sys- 
tems: Whereas, for instance, in a double-barrier scat- 
tering problem, we need to increase the energy over two 
consecutive states to re-arrive at the starting point, for 
the wire with a side branch an energy increase over one 
state only is sufficient to bring the system back to the 
same point. We also remark that Fig. is composed of 
combinations of the paths shown in Figs. |l| and |5|. 

Clearly it would be desirable to classify all the possible 
paths that are taken in the complex plane as we proceed 
from one transmission peak to another. Here we have 
emphasized only two paths, namely those shown in Fig. 
and 0. These two possibilities are a direct consequence 
of the condition Eq. ( |13|) , but of course there might exist 
scattering problems which do not obey this condition. 

To summarize this Section: We have shown that the 
addition of a particle to the system can lead to at least 



two very distinct paths of the transmission amplitude in 
the complex plane. For paths of type A (see Fig. [|) 
the parity of the Aharonov-Bohm effect in a two ter- 
minal geometry is out of phase on consecutive resonant 
peaks, whereas for a scatterer of type B (see Fig. ||) the 
Aharonov-Bohm oscillations are in phase. With respect 
to the experiment of Yacoby et al. it is an in-phase be- 
havior that is needed to explain the data. 



FIG. 6. Current amplitudes in the double-barrier model. 

We assume that the two potential barriers in this model 
are identical and that each potential barrier is symmetric 
with respect to particles incident from the left and the 
right. The current amplitudes aj, bj, Cj and dj, j = 1,2 
which are shown in Fig. pL are connected by the following 
relations, 
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FIG. 5. Schematic representation of the path of the the 
transmission amplitude as a function of the particle energy 
for the case that there is a transmission zero. The energy 
interval covered is that needed for the transition from one 
resonant peak (in the transmission probability) to the next 
resonant peak. The phase changes by about 2n. 



IV. EXAMPLES 

In this section we present the calculations which lead 
to Figs. [I], H and |3| for the transmission amplitude. In 
addition we present results for the density of states as a 
function of energy which we compare with the transmis- 
sion probability. 



A. The resonant double-barrier 

The first example is a one-dimensional double-barrier 
model, which consists of two consecutive potential bar- 
riers which scatter particles moving along the cc-axis. A 
schematic illustration is shown in Fig. o. 
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Here, f and i are the reflection and the transmission am- 
plitude of the potential barriers, respectively. The am- 
plitudes in the well are related by b\ — rc2 and c\ = rb-2 , 
where r = exp(iip) is the transmission amplitude of the 
well and ip the phase accumulated by a one time traversal 
of the well. 

We have to find the scattering matrix S which relates 
the outgoing current amplitudes 02,^2 to the incoming 
current amplitudes a\,d\. A little algebra gives 
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To proceed we parameterize the scattering matrix of the 
single barrier also in terms of angular variables (see Eq. 
(Q)). We take f = ie 10 sm.<j) and i = e 10 cos(f> with real 
angles 9 and (p. In our model the potential is uniform 
(except for the two barriers) and if is thus connected 
to the energy £ of a scattering particle via if = kl 
where I is the distance between the potential barriers 
and k = y/2mE/h is the wave vector of the particle away 
from the barriers. We assume that the quantities 9 and <f> 
are independent of the energy E. This implies in partic- 
ular that there is no particle density inside each potential 
barrier. With these specifications the transmission am- 
plitude t is given by 



t = 



el (kl+28) CQS 2 ^ 
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(18) 



The transmission amplitude is a function of energy only 
through the energy dependence of the wave vector k. 
This result is used to give the transmission amplitude in 
the complex plane as a function of energy in Fig. pi. The 
parameters chosen for Fig. [l] are 9 = 2.2 and (f> = 2.1. 
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B. Phases and density of states in the double-barrier 

model 



For comparison with the wire connected to a side 
branch, we now examine the phases in the double barrier 



scattering problem and the density of states. To find the 
Friedel phase we use Eq. (0) and find from Eqs. (O), 
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The derivative of the phase 9 {t> of the transmission am- 
plitude is found from tan# (f) = Im{£}/Re{£}, 
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Substituting Eq. (|18|) into Eq. (|20[) we can show that 
for the resonant double barrier the derivatives of the two 
phases are identical, 
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This result just restates Eq. (|fl]) for the double-barrier 
model. The phase 6 = 9 (t> = 9 {S) - tt/2 for the dou- 
ble barrier model as a function of kl (wave-vector times 
well width) is shown in Fig. Q]. The specific parameters 
arc 9 = 2.2 and <f> = 2.1. Thus the phase shows a step 
like behavior. It is nearly constant as a function of k 
and increases sharply when k, respectively the energy E, 
coincides with a resonant state. To show this, we now 
examine the density of states. 

For a perfect one-dimensional wire the density of states 
per unit length and in a small interval of wave vectors is 
dn/dk — l/(27r) for carriers moving to the right. Since 
dE/dk — hv where v is the velocity of the carrier, the 
density of states per unit length and in a small energy 
interval is v = dn/dE = l/(hv). We are interested in the 
density of states in the region between the two barriers. 
Denoting the scattering states with unit amplitude of in- 
cident carriers from the left by $i(x) and the scattering 
state of carriers with unit amplitude of carriers incident 
from the right by, $2(2;), the local density of states can 
be expressed byEil 
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The density of states in the well region (between the two 
barriers) is the integral over the local density of states 

p = ^_ l , 2 dxv(x). With an explicit calculation of the 
scattering states, we obtain 
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Here, for simplicity, we assume that the phase increment 
of a barrier traversal is zero. We have considered here 
the density of states which are obtained in terms of en- 
ergy derivatives of phases or scattering matrices (as it is 
widely done). Such derivatives do not naturally lead to 



an answer for a spatially local density of states. The den- 
sity of states in the region between the two barriers, is 
however, a local question. A rigorous procedure to obtain 
local densities is via derivatives of phases and scattering 
matrices with respect to local potentials.El The discussion 
given here (in terms of energy derivatives) is correct only 
up to WKB like corrections. The last term in Eq. ( |23| ) 
contains a factor l/(kl) and is thus small for wells which 
are much larger than a wave length. Neglecting this term 
we obtain, 
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FIG. 7. Phase of the transmission amplitude as a function 
of kl(— W2mE/h) for the double-barrier model. 

As is well known, for the case of opaque barriers (9 — > 7r/2 
and </> — * tt/2) the density of states becomes a series of 
delta functions limf^_i p = Y^^Sl 8(E—ir 2 h 2 n' 2 /(2ml 2 )) 
which coincide with the peaks of perfect transmission. 
For wells much wider than the wavelength, we find from 
Eq. (p4) that the peaks in the density of states are at 
the energies 
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Here n takes any integer values larg er than (29 — ir)/(2ir). 
Comparison of Eq. ( |l9| ) with ( |24|) leads to d9 u) /dE = 
■n p, i. e. the Friedel sum rule in the double-barrier model. 
Using Eqs. (O) and (J23) we obtain 
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a relation between the density of state p and the trans- 
mission probability \t\ 2 . This relation implies that the 
energy values £ n which determine the peaks in the trans- 
mission probability coincide with the energies E n of the 



peaks in the density of states. Using this fact and Eq. 
(E4T) we can show 
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In particular, Eq. (fL3|) applies to the resonant double- 
barriers with well widths large compared to the wave 
length. Fig. shows the transmission probability and 
the density of state as a function of kl for the double- 
barrier for the parameters 9 — 2.2 and <j> — 2.1. We 
emphasize this behavior of the double barrier, since as 
we show in the next section, a wire with a side branch, 
exhibits a transmission probability and a density of states 
which do not peak at the same energy. 
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FIG. 8. Transmission probability (the dashed line) and the 
density of states p (the solid line) in units of po = ml 2 / (nTi 2 ) 
as a function of kl(= ly/2mE/ti) for the double-barrier. 



C. Wire with a side branch 



The second example which we investigate is a perfect 
one-dimensional wire with a side branch which we also 
call the stub model. Such models have already a long his- 
tory. The conductance of a wire with a side branch was 
investigated in Refs. |34|,p5| and EJ8l Refs. |3l and |37] con- 
sidered charging effects in structures with side branches. 
In the context of the present work the phase behavior of 
the transmission amplitude has been investigated by Deo 
and Jayannavanl3i23 and Ryu and ChoEj However, nei- 
ther of these two works makes the distinction between the 
Friedel phase and the phase of the transmission ampli- 
tude. Below we investigate these two phases for a perfect 
wire to which a side branch of length V is attached. A 
schematic illustration of this system is shown in Fig. 0. 
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FIG. 9. Current amplitudes for the model of a wire with a 
side branch. 

The junction between the wir . e and the side branch is 
described by a wave splitter .E3H We consider the time- 
reversal invariant case only and consider a splitter which 
is symmetric with respect to carriers incident from the 
left and right lead. Furthermore we assume the poten- 
tial away from the junction is the same in all branches. 

The current amplitudes a'-, &'■ and c'„, j = 1,2, which are 

n 3 3 J J ] 

shown in Fig. H are connected by the following relations, 



(28) 



and b[ = r'b' 2 = e itp 'b' 2 with ip' = 2kl' + it. Here k is the 
wave vector of a particle with energy E. The constant 
■k in tp guarantees that the wave function has a node at 
the upper end of the stub. In Eq. p3), f' is the reflec- 
tion amplitude from a lead to itself, t' is the transmission 
amplitude from a lead to another lead, a is the reflection 
amplitude from the stub to itself, e is the transmission 
amplitude from a lead to the stub or from the stub to a 
lead, and r' is the transmission amplitude by which the 
particle starting from the junction returns to the junction 
through the stub. 

For the scattering matrix, relating the incident ampli- 
tudes on the wire a' 1: cK to the out-going amplitudes a' 2 , c 2 




in the wire (see Fig. |9|) we obtain 



V ?' 



1-ctt' 



1 1 
1 1 



(29) 



Here the first matrix on the right hand side arises from 
direct transmission past the side branch and direct reflec- 
tion at the wave splitter due to the side branch, whereas 
the second term of the right-hand side of is the contri- 
bution due to carriers which enter the stub and thus a 
undergo multiple scattering process. We assume that the 
scattering amplitudes f' , t' , a and e are real numbers, 
e 7^ and are independent of the energy E. These as- 
sumptions, and the fact that the scattering matrix Eq. 
(Elf) must be unitary, demands^ 



Ai+A 2 Vl-2e 2 )/2 



(30) 



i'= f-A 1 + A 2V / l-2£ 2 )/2 



o = -\ 2 \J\~2e 2 



(31) 

(32) 



where \j is 1 or —1 (j = 1,2). Depending on the 
choice of the X's four different wave splitters are ob- 
tained. The value of the coupling constant e is in the 
interval [— l/\/2, l/v2]. Using such a wave splitter leads 
to a transmission amplitude 



likl' 



t = 



-A 1 + A 2 Vl-2e 2 1 + Aie : 

2 „2ikl' 



1 - A 2 \/l - 2e 



(33) 



The path of this amplitude in the complex plane as a 
function of energy is shown in Fig. |]for the case Ai = —1. 
The path is a circle through the origin since Eq. ( p3| ) im- 
plies I t + Ai/2 | 2 = 1/4. From Eq. (||) it follows that 
the wire with the stub has zero points of the transmission 
probability at the energies 



E< 



7T 2 h 2 

2ml' 2 



1 + Ai 



(34) 



n = 1,2, •••. We re-emphasize that in contrast to the 
case of the double barrier, the origin is included in the 
path of the transmission amplitude. 



D. Phases and density of states in wire with a side 
branch 

Let us now investigate the Friedel phase and the phase 
of the transmission amplitude for the wire with a side 
branch. Using Eqs. (J7|) and (|29|), the derivative of the 
Friedel phase 9 U) with respect to the energy E in the 
wire with a side branch is given by 



d0 ( » 

dE 



dip' 
dE 



(35) 



As shown in the previous subsection there are zero points 
in the transmission amplitude as a function of the en- 
ergy in the stub model, so that abrupt phase changes of 
the transmission amplitude do occur. Thus we need a 
limiting procedure to define the phase of the transmis- 
sion amplitude. To this end we add a small perturbation 
±77, 77 > to the transmission amplitude 



i. 



t±r) 



(36) 



and evaluate the phase of t± in the limit 77 — > +0. The 
derivative of the phase 9 {t) of the transmission amplitude 
with respect to the energy E is thus given by 



QQW 

~~d¥ 



lim 



/ lmfei yl d Imfe} 
^Re{t ± }y j 8ERe{t±}' 



(37) 



Using our specific result for the transmission amplitude 
we obtain 



d0<- 



d9 ( 



+ OG 



dE 



dE 



± Ai n Y^ S ( E - E 1 -" 



(38) 



A detailed derivation of Eq. (|38j) is given in Appendix A. 
Fig. [n] shows the phase of the transmission amplitude of 
the wire with a side branch as a function of kl' . We have 
chosen the branch of the wave splitter with Ai = —1, 
A 2 = 1 and a coupling constant e 2 = 0.35 
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FIG. 10. Phase of the transmission amplitude #'*' (solid 
line) and Friedel phase 0" ' — it/2 (dashed line) as a function 
of kl' for the wire with the side branch. 

Let us next investigate the density of states. The wire 
is taken to be on the x-axis, with the splitter located at 
x = 0. The stub points along the positive y-axis. The 
splitter is described by an energy independent scattering 
matrix and can thus be viewed as point like. A scatter- 
ing state of unit amplitude expji/c^x} describing particles 
incident from the /x-th lead (/i = 1,2) gives in the side 
branch rise to a wave ^(y). Here k\ = k for a wave 
incident from the left and fc 2 = — k for a wave incident 
from the right. The local density of states v{y) in the 

side branch is given by v(y) = Y^=i I * ' niv) ? /(M, 
and the total density of states is given by the integral of 
the local density of states over the entire length of the 

stub, p — L dyv{y). We find 



V_ 
hv 



2s 



1 - ae i{ P' 



sin ip 

if' — 71 



(39) 



In the WKB limit of interest here, for a side branch much 
longer than the Fermi wavelength, the second term in 
the bracket of the right-hand side of Eq. (J3^) can be ne- 
glected. Using Eq. (32) we obtain for a long side branch 
the density of states 



/' 



As 2 



hv 



1 - VT - 2e 2 ) +4Vl-2e 2 ■ sm 2 (K 2 l') 



(40) 



Here, Kj,j = 1, 2 is defined by K = k+ (1 - \j)n/(W). 
For a stub much longer than a wavelength, the energies 
at which the density of states peaks are given by 



EL 



,2*2 



2m/' 2 



1- A 5 



(41) 



In the week coupling limit e — ► +0, using Eq. (MQ), we 

obtain lim 6 ^ +0 p' = J2 n =i^(^ ~ ®'n)- Here the right 
hand-side is the density of state of a particle confined in 
the completely isolated stub but taking into account that 
the phase change at the closed coupler is as dictated by 
A 2 = 1. The comparison of Eq. (§|) with Eq. (jl(j) leads 
to d0 u) /dE — irp'. We have thus verified the Friedel 
sum rule for the wire with a side branch. 

Let us now show that in this scattering problem (for 
the splitters with A1A2 = —1) the peaks in the trans- 
mission probability do not co-inside with the peaks in 
the density of states. Using Eq. ( p3[ ) the transmission 
probability | t | 2 can be expressed in the form, 



t? 



(l-A 1 A 2V / l-2e 2 ) cos 2 (iCJ') 



(l-AiA 2 Vl-2e 2 ) +4AiA 2 \/l-2e 2 -sin 2 (iri/') 



(42) 



Therefore the energy values £ n , n = 1, 2, • • • at which the 
transmission probability peaks are given by 



2*2 



f — 



n 2 h- 
2mV 



1-Ai 



(43) 



We see that for AiA 2 = —1 (i. e. depending on the choice 
of the splitter), these energy values are not equal to the 
energies E' n , n = 1, 2, ■ • • which determine the peaks of 
the density of states. Indeed, using Eqs. ( p0| ) and (E^) 
we obtain for the relation between density of states and 
the transmission probability 



t? - 



1 
Aid 



t' A - 



nhve 2 
IV 



(44) 



Using 



Eqs 
+1 dE 



(EO) and (Ejl) we can show that 
p' = 1. This implies that the condition Eq. 



) is fulfilled also by a wire with a side branch. The 
different behavior, depending on the sign of A1A2, is espe- 
cially apparent in the week coupling limit. In this limit, 
for AiA 2 = — 1, almost all the particles incident from the 
wire on the wave splitter pass through the junction with- 
out noticing the side branch. Conversely, in the week 
coupling limit of a wave splitter with AiA 2 = 1 almost 
all the particles incident from the wire on the wave split- 
ter are reflected at the wave splitter. Fig. [ll] shows the 
transmission probability and the density of state as a 
function of kl' for the wire with a side branch connected 
by a junction with Ai = — 1, A 2 = —1 and e 2 = 0.35. 
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FIG. 11. Transmission probability (the dashed line) 
and the density of states p (the solid line) in units of 
p' = ml' /(2irh ) for a wire with a side branch. 



E. Wire with scattering and a side branch 

The previous two models are examples which demon- 
strate two different behaviors of the transmission ampli- 
tude in the complex plane. These different behaviors are 
illustrated by Fig. [| and |5|. Clearly, both of these models 
are very particular (non-generic) and the question arises 
how the behavior exemplified by these two simple models 
shows up in more complicated structures. To examine 
this question we now consider a structure which incor- 
porates both the resonant double barrier and the side 
branch. A schematic illustration of this system is shown 
in Fig. 0. 
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FIG. 12. Current amplitudes for a model with scattering 
in the wire and with a side branch. 

We use the same potential barrier as in the resonant dou- 
ble barrier structure and the same wave splitter as we 
used for the description of the wire with the side branch. 
Again we will assume that the potential outside these 
scatterers is everywhere the same. In this model the cur- 
rent amplitudes a'L b'L c", d" and e'j, j = 1, 2 which are 
shown in Fig. [13 with the directions, are connected by 



the following relations, 



"2 



b'{ 



f t 
i f 

r 
r 

e V 



b'{ 



d'i 




(45) 
(46) 

(47) 



Furthermore, as in the model with the side branch we 
put d'l = r'd 2 '. 

For the overall scattering matrix of these system we 
obtain 



S = 



1 



V 



w 



D \ W -AitW* 



in which D, V and W arc defined by 

D = 1 - (xr' - ff (1 + Ait')t 2 
V = £>- I £ | 2 (1 - err') 
W = f*£i'(l - Ait')t. 



(48) 

(49) 
(50) 
(51) 



The transmission amplitude £ can be brought into the 
form 



tt'(l + X ie 2lkl )e 



2ikl' \„ikl 



1 + cre 



2/7,7' 



— ff'(l — Xie 2lkl ')e 



2ikl'\„2ikl ' 



(52) 



In this representation the transmission amplitude t de- 
pends on the energy E only through the wave vector 
k = \/2mE/h. As a function of energy, the path of this 
transmission amplitude in the complex plane is shown 
in Fig. [| Here, the parameters chosen for Fig. || are 
6 = 2.2, 4> = 2.1, Ai = -1, A 2 = 1, e 2 = 0.35, I = 1 
and V = 4. This figure shows clearly that a more gen- 
eral model combines the behavior of the resonant double 
barrier model and the stub model. Sequences of turns 
of the transmission amplitude which path through the 
origin are interrupted by "double turns" characteristic 
of the resonant double barrier in which the transmission 
amplitude is non-zero. From Eq. (p3) it can be noted 
that the wire with scattering and with a side branch has 
the same zero points for the transmission amplitude as 
the wire with the side branch. 



the transmission amplitudes. Except in special cases, 
these two phases are in general different- This important 
point has also been emphasized by Lee.EO The two phases 
are different if the transmission amplitude exhibits an 
energy at which it is zero. At these energies the trans- 
mission amplitude phase and the Friedel phase acquire 
an additional difference given by ±7r. If the transmission 
amplitude exhibits no zero point between resonances the 
variation of the phase from one resonant peak to another 
is close to it. If a zero point exists the phase change 
between consecutive resonance peaks of the transmission 
probability is close to 2n. This difference is shown clearly 
in the paths of the transmission amplitude in the complex 
plane. 



For a sufficiently general model, we expect sequences 
of resonant peaks which are in phase (the phase increases 
by 27r) interrupted by peaks which are out of phase (the 
phase increases only by 7r as we go from one peak to 
the next). In terms of the parity of the Aharonov-Bohm 
effect, these implies sequences of peaks over which the 
parity is conserved. These sequences are interrupted by 
transitions which generate a flip in the parity. It is clearly 
desirable to investigate now a number of statistical ques- 
tions: For example for a fully chaotic quantum dot one 
would like to find the ensemble averaged density of zero's 
and compare this with the ensemble averaged density of 
states. Furthermore, one would like to know if such a 
cavity exhibits correlations in the occurrence of zero's 
(long sequences of zero's interrupted by a flip in the par- 
ity), etc. Since the distribution of eigenvalues is known, 
random matrix theory likely gives an answer to these 
questions. 



We add a remark on Fano resonances: Fano resonances 
arise due to the coupling of discrete states with contin- 
uous states. Such resonances also exhibit transmission 
zero'sOCj Consequently, for such resonances the Friedel 
phase also does not coincide with the phase changes of the 
transmission amplitudes. The wire with a side branch in- 
vestigated here also couples a set of discrete states with a 
continuum and thus provides for interfering transmission 
paths. However, the resonances in this case are not of the 
Fano type (as shown in Fig. [by) but rather Breit-Wigner 
resonances in the reflection probability. 



V. CONCLUSIONS 

In this work we have discussed the transmission am- 
plitude as function of energy in the complex plane for 
scattering systems (without a magnetic field) connected 
to two single channel leads. We emphasize that there are 
two phases which are of importance, namely the phase 
which appears in the Friedel sum rule and the phase of 



It is very interesting to investigate the behavior of the 
transmission amplitude in the complex plane for varia- 
tions in parameters other than the energy. For instance 
we can ask about the path of the transmission amplitude 
if the AB flux increases by a flux quantum in a multiply 
connected geometry. The additional questions raised in 
this section clearly demonstrate that the investigation of 
the path of the transmission amplitude in the complex 
plane is an interesting avenue of future research. 
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APPENDIX A: RELATION BETWEEN PHASES 
FOR THE WIRE WITH A SIDE BRANCH 

In this Appendix we outline the derivation of Eq. (p3) 
which relates the Friedel phase and transmission ampli- 
tude phase for the wire with a side branch. First, using 
Eqs. © and Q we find 



i± 



-Ai — lt 1 — Ai e 



%ip 



1 - cr e % f 



T±V- 



(AI) 



The transmission amplitude t is dependent on the en- 
ergy only through the quantity ip', so that we obtain 



d8 {t) 



lim 



<V 9Arg{i±} 



dE n-*+o dE dip' 



(A2) 



Using Eq. (AI) the derivative dArg{t±}/dip' is found to 
be 

dArgfe} 
dip' 

I \Rc{i ± }J f %'Rc{t ± } 



i 



\-G l 



2 1 + a 2 -2a cos <p' 



T2fi(^cp')v + h(<J,\up') 
2f 2 (a,ip')v 2 + (lT2X 1 r 1 )f 3 (a,X 1 ,ip') 
1 1-a 2 



2 1 + cr 2 - 2cr COS <£' 
1 



-FM T 2Al7 7 , AlCT.^' + i^-TT 



(A3) 



where the functions fi(x, y), fz(x,y) and f^{x,y,z) are 
defined by 



fi{x,y) = 2x- (l + x 2 )cosy 

f 2 (x, y) = 1 + x 2 - 2x cos y 

h{x,y,z) = (l + xy) 2 (l -ycosz) 

and the function F(x, y, z) is defined by 

F(x,y,z) 

s (l-x 2 )(l-y 2 ) 

(1 -x) 2 (l -y) 2 + 2(1 + xy)(a; + y)(l + cosz)' 



(A4) 
(A5) 
(A6) 



(A7) 



F(x, y, z) has the following properties. First it follows 
that 



lim F{x,y,z) 

in z = (2n + 1)tt 

+e(l — y 2 ) x oo inz^ (2n + l)7r, 

n = 0, ±1, ±2, • • •. Second it follows that 

I-X+2TT 



(A8) 



dzF(x,y,: 



— <b du>- 



(l-x 2 )(l-y 2 ) 



i Jc {(I + xy)uj + x + y} {(x + y)uj + 1 + xy} 

(A9) 



-27T in ^ | x + y \<\ 1 + xy \ 
-2n in \ x + y \>\ 1 + xy \ ^ 



where A is a real number and ix> = exp(zz). C is the path 
running counterclockwise on the cir cle w hose cen ter is 
the origin and the radius is 1. Eqs. (A8) and ( |A9| ) lead 
to 



lim F(x,y,z) 

x— *1±0 



+oc 



= +27re(l-y 2 ) ^ 6(z - (2n + l)vr). (A10) 

n— — oo 

Using Eq. ( J34| ) and the inequality _E > 0, we obtain 

(AH) 



From Eqs. JA2| ), ( |A3| ), ( |A10| ) and ( |ATl| ) we derive 



86^ _ dp' 
~d~E ~ ~dE 



1-cft' 



+ 00 



±Xnr^2S(E-E^). (A12) 



Using Eqs. (|35D and (|A12|) we obtain Eq 
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